Quench dynamics of the interacting Bose gas in one dimension 
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We obtain an exact expression for the time evolution of the interacting Bose gas following a quench 
from a generic initial state using the Yudson representation for integrable systems. We study the 
time evolution of the density and noise correlation for a small number of bosons and their asymptotic 
behavior for any number. We show that for any value of the coupling, as long as it is repulsive, the 
system asymptotes towards a strongly repulsive gas, while for any value of an attractive coupling 
the long time behavior is dominated by the maximal bound state. This occurs independently of the 
initial state and can be viewed as an emerging "dynamic universality" . 



The interest in quantum systems out of equilibrium 
has been re-ignited with the advent of ultra-cooled atoms 
in laser traps or lithographic nano-devices, which allow 
finely controlled experiments to explore various aspects 
of non-equilibrium dynamics - quenching or non-linear 
transport being examples In parallel, advances 

in theoretical computational techniques now allow us to 
calculate the time evolution of various observables of in- 
teracting many body systems and study the underlying 
physics Many of the systems studied are governed, 

to a very close approximation, by integrable Hamiltoni- 
ans 0, E| , opening the way to understanding their quan- 
tum dynamics in great detail. 

We shall be concerned here with non-equilibrium 
quench dynamics - preparing the system in some initial 
state |3>o) and following its evolution under the influence 
of a Hamiltonian H applied to it suddenly at time t = 0. 
For example, a Bose gas initially in a Mott state in a deep 
periodic trap, is time evolved under a supcrfluid Hamil- 
tonian with the trap depth reduced, allowing the bosons 
to easily hop from site to site. The evolution can be for- 
mally obtained by expanding the initial state in terms of 
a complete set of orthonormal energy eigenstates {A}, 
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For an exactly solvable model, these eigenstates |A) are 
given by the Bethe Ansatz. This approach has been effec- 
tive in obtaining the full spectrum (and thermodynam- 
ics) of a number of important models (3-Hlj. However, 
for time evolution, in addition to the spectrum, one also 
needs to compute the overlap of the eigenstates with the 
initial state and carry out the summation - difficult tasks 
due to the complexity of the Bethe Ansatz eigenstates. 
In spite of this, some progress has been made 12Ml4|. 

In this letter, wc address this question generalizing a 
"contour integral" approach introduced by V. Yudson for 
the purpose of studying superradiance effects in the con- 



text of the Dicke model [15|, [16[ . The contour represen 



tation does away with both the above steps, and replaces 
the summation over eigenvalues with an integral along 
contours in the complex plane, appropriately chosen to 



faithfully represent the initial state, capturing the rel- 
evant overlaps via residues of poles in the two particle 
S-matrix of the Bethe Ansatz eigenstatc. 

We show that the framework has general applicabil- 
ity and use it to understand the quench dynamics of an 
interacting gas of bosons, with short range attractive or 
repulsive interactions, 

H= [ db\x)db{x)+ctf{x)b{x)tf{x)b{x). (2) 



We shall obtain the time evolution of an arbitrary initial 
state, and calculate the evolution of the density and noise 
correlations. 

The model ([2]) was solved by Licb and Liniger Q. Its 
N- particle eigenstates, labeled by the momenta {A}, are: 
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where S y is a symmetrizcr. 



The factor, ZfAz) 



lC S z-ic' Vj > mcor P ora tes the S-matrix, Sij (A, — Xj) = 
that describes the scattering of two bosons with 
The corresponding eigenenergies are 



Xi — Xj -\-ic 
Xi — Xj 



momenta A; , Xj 



m E\ = Aj with the momenta being real- valued for 



repulsion (c > 0), or complex conjugate pairs (signifying 
bound states) for attraction. One usually proceeds by 
imposing periodic boundary conditions to determine the 
values of the {A}'s and subsequently studies the infinite 
volume limit. Instead, we shall work directly in the infi- 
nite volume limit with the momenta {A} unconstrained, 
allowing us to integrate over them instead of summing 
over discrete values. 

A given state |$ ) = J x $ (xi ■■■x N ) I\f=i b H x j)\ Q ) 
($o symmetric) can be directly time evolved if it is ex- 
pressed in terms of the eigenstates, eq. It can be 
shown that, (we denote 6{x\ > ■■■ > xn) = 
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where the momenta integration contours are determined 
by the interaction. In the repulsive case, all {Xj} have to 
be integrated along the real line, while in the attractive 
case, the Xj have to be integrated along lines parallel to 
the real axis that are separated in the imaginary direc- 
tion by a little more than |c|. This separation captures 
the bound states that appear in the spectrum of this 
model. The proof of the statement involves studying or- 
dered initial states, \x) = 9{x) Y[j b^(xj)\0) which, when 
represented in the form above, have zero residues for the 
chosen integration contours. It can then be shown that 
each Xj integration reduces to S(yj — Xj) proving the va- 
lidity of (@|. For details of the formalism as applied to 
the interacting Bose gas, please see 
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In the rest of this article, we discuss the quenching of a 
system of bosons initially in a Mott state (a deep periodic 
trap) or in a condensate (a parabolic trap) with the trap 
removed at t = and the system evolving subsequently 
under the action of the interacting Hamiltonian ©. We 
shall be interested in the effects of interaction, attractive 
and repulsive, on the evolution. 

We first discuss the case of bosons trapped in a periodic 



potential |<I> 
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If we assume that the wave functions of neighboring 

a 2 

bosons do not overlap significantly, i.e., e~~^ <C 1, then 
an ordering of the initial particles (required in the Yud- 
son representation), with Xj = (l—j)a, is induced by the 
non-overlapping support. To illustrate this, we find the 
exact two particle finite time wave function: 
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where a(t) = 2ct — i(yi — x\) — i(y 2 — x 2 ). In the attrac- 
tive case with the same initial conditions we find the same 
result with erfc (the complementary error function) re- 
placing erf. This induces a significant change in behavior, 
as can be seen from the evolution of the density. Fig. Q] 
shows the evolution of the density measured at x = 
for free, attractive and repulsive bosons. While the re- 
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pulsive, and non-interacting cases are nearly identical the 
attractive case shows oscillations with period T ~ , due 
to the competition between diffusion of the bosons and 
the action of the attractive potential. This can be seen 
in the time evolution of the wave- function for two attrac- 
tive bosons, (Fig. [2]). As time evolves, the initial peaks 
coalesce and increase in height. At longer times, this cen- 
tral peak "breathes" as it diffuses. The oscillations arise 
due to the formation of a bound state that appears as a 
contribution from a pole of the S-matrix [181 ]. These ef- 
fects are measurable in time-of-flight experiments in both 
cases. We expect the results to be qualitatively similar 
for any number of particles. 




FIG. 2: {p(x, t)) vs. x of two attractive bosons after a quench 
from |</>iatt) plotted for different times. The two separate 
peaks coalesce and form a bound state as t — > oo. The key 
shows the times (in units of a 2 ). 

Exact integrations beyond two particles are difficult. 
However, the contour approach allows us to extract the 
asymptotic behavior of the evolved wavefunctions. Con- 
sider repulsive bosons and carry out the scaling A — > Xs/i. 
Then, Zf- (Xi - Xj) -t sgn^ - yj )[l - ^(A,- - A,-) + . . .], 
yielding to leading order, 
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c'(y) being fcrmionic creation operators replacing 
the "fermionized" bosonic operators, Yij c illj) — 

Ui<j sgn(y* - yj)b\yj). H$ = J x dc){x)dc(x) is the free 
fermionic Hamiltonian. A y is an anti-symmetrizer act- 
ing on y. Thus, the repulsive Bose gas for any value of 
c > 0, is governed in the long time by the c = oo hard 
core boson limit (or its fermionic equivalent) 
and the system equilibrates with an asymptotic momen- 
tum distribution, nk = ($o|cJ,Cfe|$o), determined by the 
antisymmetric wavefunction <&o{y) = •A y 9(y)(&o{y) and 
the total energy, E<s, = ($o|-ff|*o)- 

The corrections to the asymptotics, can be obtained 
using the stationary phase approximation to carry out 
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the A integrations. We have for the repulsive case, 
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where £ y = Mr and we drop ^ since the initial condi- 
tions have finite spatial extent. Asymptotically, the evo- 
lution of the density is not instructive, so we compute 
the iV-boson density-density (noise) correlation function, 
P2(z,z') = (p(z)p(z')). At long times, 



p 2 (z,z';t) -> p 2 (£z,£z>) = 
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with g zz , = l-2cV^cr5(^-C 2 '-ic)[e( c + l «^ 2<T2 erfc((c + 

i6)o-) + e( c -^') 2<T2 crfc((c-^ 2 0cr)] and^ = z/2t. This 
quantity was also studied by Lamacraft [2l| for the re- 
pulsive model. However, our result differs from his. 

An interesting way to observe the evolution of repul- 
sive bosons to free fermions is the Hanbury-Brown Twiss 
(HBT) effect. Consider the normalized spatial noise cor- 
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P2(z,z') 



p(z)p(z') 



1 



C2(z,z'). In the non-interacting 
case S(£) = 1 and g zz i = 1 and we recover the HBT 
result for N = 2 Q, C 2 °(£ 2 ,£ z ') = ±cos(a(£ z - &,))• 
At any finite c, we can see a sharp fcrmionic dip appear 
(anti-correlation) that gets broader with increasing c as 
shown in figure [3] For higher particle number, as shown 
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FIG. 3: Normalized noise correlation function C2(£, — £). 
Fermionic correlations develop on a time scale r ~ c -2 , so 
that for any c we get a sharp fermionic peak near £ = 0, i.e., 
at large time. The key shows values of ca. 

in fig. 31 we see "interference fringes" that get narrower 
and more numerous with an increase in number of parti- 
cles. However, the asymptotic fermionic character is still 
visible. 

The scaling argument fails for the attractive g as as 
the As are integrated on parallel trajectories, and as t 
increases, their separation also grows. The asymptotic 
behavior of attractive bosons, as we shall see (and already 
observed earlier), is dominated by a maximally bound 
state. Since the contours of integration spread out in the 




FIG. 4: Normalized noise correlation function — £) for 

five repulsive bosons released from a Mott-like state for ca — 
2. 




FIG. 5: Contribution from stationary phase and pole at large 
time in the attractive model. The blue curve represents the 
shifted contour. 



imaginary direction, in addition to the stationary phase 
contributions at large time, we also need to take into 
account the contributions from the poles. Figure [5] shows 
an example of how this works. Taking into account all 
the pole contributions, we get a sum over several terms 
given by, 
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with Zij(-ic) = c^/t6(yj > j/j). Note that while the 
asymptotic dynamics of the repulsive model is given 
purely in terms of the "lightcone" variables £ 3 = ||, 
this is not the case in the attractive model where observ- 
ables acquire explicit time dependence. Also, important 
to note is that the bound state contributions which ap- 
pear when one or more of the S-matrix poles are picked 
up, are higher order in t as compared to the non-bound 
states. Asymptotically therefore, the larger bound states 
dominate. In fig. |6] (obtained by numerically integrating 
the expression for the noise correlation) we see interfer- 
ence fringes similar to the repulsive case with the central 
peak increasing positively indicating the dominance of 
the maximally bound state. 

We now consider the evolution of the Bose gas af- 
ter a quench from an initial state where all the bosons 
are in a ground state of a harmonic trap, |$ C ond) = 
f FT — ^re-^/^bUxiMO). In order to use the Yud- 

son representation, we have to rewrite the initial state 
as a sum over all different orderings of the coordinates 
using the symmetry of the initial wavefunction $o- For 
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FIG. 6: Ca(£, — £) for three particles in the attractive case 
plotted for three different times. At larger times, the cor- 



relations away from zero fall off. 
magenta and yellow respectively. 
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20, 40, 60 for blue, 



the repulsive model, the stationary phase contribution is 
all that appears: 
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For the attractive shown in eq. (|5|), the bound 

states dominate, with the larger bound states (more pole 
contributions) dominating at larger time. This holds in- 
dependently of the initial conditions. 

Figure [7] shows the noise correlation for two and three 
repulsive bosons starting from a condensate. We see the 
characteristic fermionic dip develop. The plots for three 
attractive bosons are shown in fig. [5] The oscillations 
arising from the interference of particles separated spa- 
tially (HBT) do not appear. The attractive case however 
does show the oscillations near the central peak that are 
visible in the case when we start from a Mott insulator. 
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FIG. 7: C2(£, — £) for two (blue) and three repulsive bosons 
starting from a condensate. Unlike the attractive case, there 
is no explicit time dependence asymptotically, ca = 3 

We conclude that although the details of the time evo- 
lution depend on the initial state we quench from, the 
asymptotics show universal features given in terms of 
a "dynamic RG picture" with the asymptotic dynamics 



FIG. 8: C2(£, — £) for three attractive bosons starting from 
a condensate. Note that the side peak structure found in 
fig. [6] is missing due to the initial condition. We show the 
evolution at three times. As time increases, the oscillations 
near the central peak die out. Times from top to bottom 
tc 2 = 20, 40, 60 



being controlled by c = ±oo Hamiltonians for the repul- 
sive and attractive interactions, respectively, with time 
playing the role of the successively reduced cut-off in an 
RG procedure, projecting out the effective low energy 
physics. This picture suggests "basins of attraction" of 
perturbed Hamiltonians around the Lieb-Liniger Hamil- 
tonian (e.g., a short range potential replacing the delta 
function), whose long time evolutions are given by the 
c = ±oo limit. We therefore expect our results, beyond 
their theoretical significance, to also provide experimen- 
tal predictions. 

In summary, we have used the contour integral rep- 
resentation for both the attractive and repulsive Lieb- 
Liniger model and shown that we can extract the asymp- 
totic wave function and observablcs analytically. The 
representation overcomes some of the major difficulties 
involved in using the Bethe-Ansatz to study the dynam- 
ics of some integrable systems. The results indicate the 
emergence of universal dynamic features. The applica- 
tion to other models is under consideration. 
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